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ABSTRACT 

We present a new set of analytic models for the expansion of HII regions powered by 
UV photoionisation from massive stars and compare them to a new suite of radiative 
magnetohydrodynamic simulations of self-gravitating molecular clouds. To perform 
these simulations we use RAMSES-RT, a Eulerian adaptive mesh magnetohydrody¬ 
namics code with radiative transfer of UV photons. In parallel to the simulations, 
we develop analytic models that describe the radial evolution of the HII region in a 
range of density and velocity fields. We give a radius at which the ionisation front 
should stop expanding (“stall”). If this radius is smaller than the distance to the edge 
of the cloud, the HII region will be trapped by the cloud. This effect is more severe 
in collapsing clouds than in virialised clouds, since the density increases dramatically 
over time, with much larger photon emission rates needed for the HII region to escape 
the cloud. We also measure the response of Jeans unstable gas to the HII regions to 
predict the impact of UV radiation on star formation in the cloud. We find that the 
mass in unstable gas can be explained by a model in which the clouds are evaporated 
by UV photons, suggesting that the net feedback on star formation should be negative. 

Key words: (stars:) massive, (ISM:) H ii regions < Interstellar Medium (ISM), Neb¬ 
ulae, ISM: clouds < Interstellar Medium (ISM), Nebulae, stars: formation < Stars, 
methods: numerical < Astronomical instrumentation, methods, and techniques, meth¬ 
ods: analytical < Astronomical instrumentation, methods, and techniques 


1 INTRODUCTION 

Massive stars form in dense gas clouds, where they heat 
the gas via ultra-violet (UV) photoionisation. The pres¬ 
sure difference created between these HII regions (where 
hydrogen is ionised from HI to HII) and the external 
medium causes an expansion wave that can disrupt the 
cloud in which the star formed and prevent (or en¬ 
hance, via shock compressio n) furt h er sta r formation. Pre ¬ 
vious theoretical papers bv iKa hn lfl,954 j): fSpitzei) ( 1978t); 


I Whitworth dl9T: 


( 19971) : iMatzne: 




Franco et al. 1 199(J) : I Williams &; McKee 
20021 ) have successfully described the 


broad features of such regions using analytical arguments. 
These models have been successfully reproduced in hydro- 
dynamic simulatio n s of varying physical c o mplexity by, e.g., 
IDale et al] ll2005l): iMellema et~aL (120061): iKrumholz et al] 


(l2007h: lArthur et al.l (l201ll) : IWalch et al.l ll20ll l2013h : 
IDale et alJ (120141). Applying these models to observed HII 


regions as m 


Tremblin et al.l (120141 ); iDidelon et al.l (120151 ) 
can help further our understanding of these systems, al¬ 


though there exists some degeneracy in the models that 
makes accurate comparisons difficult. 

More recen tly, iHosokawa fe Inutsukal (l2006ll and 
(j2012l ) have modified the equations given by 


Raga et al 


Spitzed ( 19781 ) for the thermal expansion of HII regions to 


include the shell inertia and external pressure te r ms re ¬ 
spect i vely. L i mits gi ven bvlOarcia-Segura fc Francol l|l99fil ): 
iKetd d2002L 120031) : IDale et all (1201211 put constraints on 
the expansion of HII regions due to pressure and accretion 
flows, leading to ultracompact HII regions that cannot 
escape the core or cloud in which they are formed. If 
the ram pressure and thermal pressure outside the region 
is equal to or larger than the thermal pressure in the 
photoionised gas, the front will be trapped. 


The behaviour 
gions is important 


of ionisation fronts and 
because of their role in 

the clouds in which they are formed . Meffema^et_alJ 


HII re¬ 
shaping 


ll200fil): lO ritschn eder et alJ ll2009l ): I Tremblin et alJ I 2012l l: 


IWalch et all ( 20121) show in numerical simulations that HII 





























































































2 S. Geen 


regions can add a significant amount of momentum to 
the Interstellar Med ium ( IS M). Further, Roger s fc Pittardl 
J2013l l: lOeen et al.1 (|2015l f : IWalch fc Naabl (I2015l f demon- 
strate that pre-supernova stellar feedback by stellar winds 
and photoionisation allows supernova explosions to propa¬ 
gate much more efficiently than if they had explod ed in¬ 
side the dense cloud (e.g. Ilffrig fc Hennebefia l2015l l. Star 
formation can also be qu enched by UV photoion i sation feed¬ 
back, a s dem o nstrated in lDale fc Bonnelll (1201 lfl : [Dale et ahl 
(120121. l2013ll . IWalch et al.l d2012h note that a small in¬ 
crease in star formation rate can be achieved as the shock 
front around an expanding HII region compresses dense gas 
clumps in a molecular cloud, but the net effect of UV radia¬ 
tion from stars on star formation in a given cloud is expected 
to be negative. 

Our goal in this paper is to use a combination of ana¬ 
lytic and numerical work to introduce a new analytic model 
for the evolution of HII regions that describe the interplay 
between self-gravity in clouds and UV photoionisation. We 
focus on the radius of the HII region with time, though this 
quantity is extended successfully to derivations of the mo¬ 
mentum and mass of the HII region. Our model is based on 
previous analytic work but is extended to describe the pho¬ 
toionisation of turbulent molecular clouds. It considers the 
balance between the thermal pressure from UV photoionisa¬ 
tion and ram pressure from turbulent or infalling gas in the 
external medium. A complete description of the dynamics 
of self-gravitating, magnetised clouds is beyond the scope 
of this paper. Instead, we adopt simple prescriptions based 
a power law density field with a velocity field equal to the 
escape velocity at each radius. We also provide models that 
sample directly the spherically-symmetric density and ve¬ 
locity profiles to compare to the simple power law model. 
To test these models, we introduce a new suite of 3D simu¬ 
lations of UV photoionisation inside turbulent, magnetised 
self-gravitating clouds with sources of ionising radiation rep¬ 
resenting massive stars. In Section [2] we present our models 
that will be applied to our simulations. We also (re-)derive 
equations for the limits to expansion that lead to ultracom¬ 
pact HII regions. In the Section [3] we introduce our suite 
of numerical simulations. Section [5] introduces the simula¬ 
tion results and addresses qualitatively the response of the 
cloud to UV photons and magnetic fields. In Section [6] we 
compare our simulations to our analytic models assuming no 
external velocity field. In Section [3 we include an external 
velocity field in our models and discuss at what point either 
model should be used. In Section[8]we consider the expected 
response of star formation in the cloud to UV photoionisa¬ 
tion. Finally we discuss the consequences of our results and 
scope for improvement in future work. 


2 EXPANSION OF AN HII REGION 
2.1 Basic Principles 

The goal of this project is to better understand the evolu¬ 
tion of HII regions in complex astrophysical environments. 
In order to interpret our results we require a theoretical 
understanding of the behaviour of HII regions in simpler 
environments that we can compare to our simulation re¬ 
sults to determine at what point our results diverge from 


existing analytic work. The simulation code used in this 
paper, RAMSES-RT, has been tested against othe r radia¬ 
tive transfer codes as well as analyt ic arguments bv ISoitzerl 
(1 1978tl ; iHosokawa fc Inutsukal (l2006tl in an idealised environ¬ 
ment and compares favourab l y to both (see the Starbench 
project, presented in Bisbas et alj[2pi5). 

In this section we review our theoretical understanding 
of the evolution of ionisation fronts and derive expressions 
that will be used to model our results. An HII region is a 
volume of hydrogen gas that has been ionised, in this case by 
ultraviolet (UV) photons. Photoheated hydrogen has a tem¬ 
perature of approximately 10 4 K. The precise temperature 
is regulated by t he de nsity and metallicity of the medium 
(see, e.g. iDraina 120111. In this work we ignore the role of 
molecular hydrogen. Its dissociation energy of 4.5 eV means 
that there will be an extra energy cost in unbinding each H 2 
molecule. In practice, however, our photons are normally 
much more energetic than the ionisation energy of hydro¬ 
gen (13.6 eV), with the extra energy being lost to radiative 
cooling, and as such we do not expect accounting for this to 
change our results significantly. 

The temperature difference between the photoionised 
and neutral hydrogen creates a pressure difference across 
the ionisation front, that triggers an expansion wave into 
the surrounding medium. The expansion of an HII region 
around a s ource o f ionising UV photons is modelled in two 
phases fsee lKahnlll954l f. In the first phase, the front expands 
rapidly but hydrostatically to an equilibrium radius at which 
the emission rate of photons from the star is balanced by 
the recombination of photoionised hydrogen inside the ioni¬ 
sation front. This is called the Stromgren radius, denoted by 
r s . This phase is complete on the order of the recombination 
time in the medium, which at molecular cloud densities is 
on the order of kyrs. Since this is much shorter than the dy¬ 
namical time of our system for our models we do not treat 
it in our models. Instead we assume that our front starts at 
the Stromgren radius, which is typically small compared to 
the radius of the HII region at later times. 

The total number of recombination events per unit time 
in a uniform, ionised sphere of radius r s dominated by hy¬ 
drogen is given by 47r/3r s n e noaB, where n e is the electron 
number density, no is the initial hydrogen number density 
and a® is the recombination rate in the photoionised gas. 
For a fully ionised gas, n e = no- Hence we can equate this 
to the ionising photon emission rate S * and write 


( 3 S t 
\47T n%a B ) 


( 1 ) 


However, because this gas is at a higher temperature than 
the surrounding cloud material, the pressure inside the 
ionised gas is higher tha n the pr e ssure outside and so the 
ionisation front expands. ISnitzerj (. 197ft ) makes the simpli¬ 
fying assumption tha t the ionisation front is c oupled to the 
shock front (see also iDvson fc William^ Il980l l. The ionisa¬ 
tion front has a sound speed a , which is set by the radiative 
cooling equilibrium. In our simulations c; = 12.5 km/s. This 
gives a c onstant temperature in the ionised gas of 8400 K 
(based on lTremblin et aI1l2014l . for an HII region at the same 
distance from the galact ic centre as ou r Sun). 

If we assume, as in ISpitzeil ( 1978 1. that the ram pres¬ 
sure on the expanding front (left hand side) is equal to the 
thermal pressure inside the ionised gas (right hand side): 
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n ext fi 2 = riiCi, (2) 

where n ex t is the hydrogen number density in the external 
medium just outside the shock radius, n is the radius of the 
ionisation front and rq is the hydrogen number density in 
the ionised gas. 

As the ionisation front expands, the gas must remain 
in photoionisation equilibrium. Similarly to the arguments 
made for Equation [T] we can write 


n 


3 S* 
4-7T nfa b 


i 

3 


(3) 


where n (t = 0) = r a . In the next subsection we solve these 
two equations for a power law density profile. 


2.2 Expansion in a Power Law Profile 

ISpitzeij dl978h assumes a uniform density outside the HII 
region. iFranco et al. 111993) instead use a power law density 
field with a flat core. We use a similar approach, but make 
the simplifying assumption that the core radius is negligi¬ 
ble, since this does not significantly affect our results. Our 
density field is thus given by 


next (r) = n 0 — 
\ro 


(4) 


where r is the distance from the centre of the cloud, and 
ro and w are parameters fit to the cloud profile in our 
simulations, with no the density at the position of the 
source, where r B is calculated. Fitting this power law to 
observed clouds is a non-trivial process, but it has been 
erformed by, e.g. lYun fe Clemensl (1 199lT) : iBac mann et al.l 
200(f) : iNielbock et al.l ( 2012)4 : Didelon et al. ( 20151) . There 
are various difficulties in turning an observed 2D column 
density map into a full 3D density profile. A spherically 
symmetric approximation is easier to measure and calcu¬ 
late, however, than a full 3D density map. We discuss the 
validity of the assumption of spherical symmetry later in the 
paper. 

We substitute n ex t in Equation [2] for its value in Equa¬ 
tion [4] We then substitute for m in Equation [3] and solve 
the differential equation in r\ to obtain 


n(t) 


( 1 (i ct Y 

4n a B ) V n o r o ) \^p C J 


(5) 


where ip = 4/(7 — 2w). The initial radius of the ionisation 
front is r s , though this is typically much smaller than the 
radius of the fr ont an d t hus ca n be neglected when calculat¬ 
ing n (see also lMatznedl2003) . For a ho mogeneous medium 
(w = 0), we arrive at the lSpitz~ei1 (ll97Slf solution (assuming 
r s to be negligible^_ 

Similarly to lMatzned (|2002i '). we can calculate the prop¬ 
erties of the HII region based on Equation [d] The mass in 
ionised gas is equal to the density in ionised gas m multi¬ 
plied by the volume of a sphere of radius n. We can use 
Equation [3] to calculate this as 


M ion 


/ 4-7T St \ 2 m-H § 

v~ 3 ~^/ 


( 6 ) 


where m,H is the mass of one hydrogen atom and A'(= 0.76) 
is the hydrogen mass fraction in the gas. 

Note that this value is typically much smaller than the 
mass displaced by the ionisation front, since Equation [3] 
evolves such that n\ oc rF 3 ^ 2 , so as the ionisation front ex¬ 
pands n\ « n ex t, provided that the power law index w is 
smaller than 3/2. This means that the mass of the shell is 
approximately the same as the mass of displaced gas, which 
can be obtained by integrating Equation [I] The momentum 
of the shell is then the rate of change of n multiplied by the 
mass of displaced gas: 


Mi' 


■ hi mu ip 4—ti 

o- r o no -r^T r i 

3 — w X t 


(7) 


In subsequent sections we compare our simulation re¬ 
sults to Equations[S]to[71 and determine at what point these 
expressions break down. In the following subsection we in¬ 
troduce theoretical arguments for why and where this should 
occur. 


2.3 Constraints on Expansion 

The above expressions assume that the only force oppos¬ 
ing the expansion of the HII region is ram pressure from the 
static neutral mass as the shock moves through it. In reality, 
there are a number of proc eses that are able to resist th e ion ¬ 
isation fron t ’s exp ansion. iGarcia-Segura fe France] (|l996i ): 
iRaga et al.l (120121) argue that the pressure in the exter- 
nal medium can cau se the ionisation front to stall, with 
iTremblin et al.l (120141 ) extending this model to the turbu¬ 
lence in the IS M at the late stages of the front’s expansion. 
As well as this. Il<eu] (|2002l ) gives some theoretical limits at 
which accretion m ust block the expa nsion of HII regions in 
stellar cores, while [Dale et al.l (l2012f) states that for molec¬ 
ular clouds if the escape velocity at the Stromgren radius 
exceeds the sound speed in the photoionised gas, the front 
remainds trapped. 

We now adapt the model given in the previous section 
to determine the point at which it breaks down in the pres¬ 
ence of an external gas infall, thermal pressure or turbu¬ 
lence. Equation [5] can be modified to include terms in tur- 
bu lence and infall v elocity (derived in Appendix [A] based 
on IRaga et al.ll2012l) : 


1 dn 

Ci dt 


where 


F{r,t) = 


= F(r,t) — 

^ext 

1 

c? 

F{r, t ) 

/ rn 

(r s 


V ^ext 

\ r i. 

) V 


+ 


next(r, t) 


( n 0 


\n e xt(r,t) 


1/2 


( 8 ) 

(9) 


Cext is a term including the sound speed and turbulent mo¬ 
tions in the gas just outside the shock radius and n ex t is 
the velocity of the gas just outside the shock radius normal 
to the shock surface (assumed in ID to be radial from the 
source position). 

In the limit where c ex t —» 0 and v ext —> 0, we arrive 
at the solution in Equation [5] For a general case where 
Cext and Vext are non-negligible, it is not possible to solve 
Equation [H] analytically. However, we can make arguments 
concerning the limit at which f\ = 0. This is the radius 
at which the ionisation front will stall and be unable to 
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expand further. In this limit for infall-dominated flows, 
F(r,t) ~ Vext J r ’ i) . Similarly, for a cloud dominated by ther¬ 
mal pressure, F(r,t) ~ iTremblin et al.l (1201 4'! use the 
approximation c ex t = l/3cr 2 where cr is the velocity disper¬ 
sion of the turbulence in the external medium. If the veloc¬ 
ity of the gas grows over time, the ionisation front can even 
shrink as the HII region is crushed by accretion. 

For a quasi-static power law density profile as given in 
Equation [4] where n = r a at t = 0, we can write 


F(r,t) = 




( 10 ) 


The radius at which the ionisation front stalls (n —> 0) is 
found by substituting Equation [TO] into Equation[5]for a flow 
dominated by infall, with an identical expression in c ex t if 
the sound speed in the neutral gas dominates: 


T’stall — 7*0 



(ii) 


where r s t a ii is the value of n at which n = 0. If w ~ 3/2, 
then the front cannot expand beyond ro regardless of radius 
if the escape velocity exceeds the sound speed. 
Alternatively, we can write 


7* stall 


( 3 

S*\ 

/ m H cf \ 

V 4-7T 

a® ) 

V X Pr, m ) 


( 12 ) 


where P ram is the ram pressure on the shock from an external 
velocity field (^n ex tu 2 xt ), turbulence or thermal pressure 
(^^rienctCext). Note that P ram and u ex t can depend on radius 
(as we discuss in the following section), in which case the 
dependence of r stall on Ci and S* will change accordingly. 


2.4 Stalled Expansion in an Accreting or 
Virialised Cloud 

The infall velocity used to calculate these limits depends 
on the dynamics of the cloud. A complete tr eatment of thi s 
subject is beyond the scope of this paper. iHunteil d 19621 '! 
gives a simple spherical collapse model in which a cloud 
accretes onto the centre at free-fall speeds. In our work, 
the cloud is supported by turbulence, as well as thermal 
and magnetic pressure. For a cloud in virial equilibrium, 
where iv 2 xt = | XX- (G is the gravitational constant, M is 
the mass of the cloud and R is the radius of the cloud), 
we can estimate the velocity of the gas inside the cloud to 

be approximately ^/|" the escape velocity, though this is a 
crude approximation to the actual velocity structure, which 
is continuously evolving. In addition these motions will be 
random, rather than directed solely towards the centre of 
the cloud. One source of infall in the cloud is from relax¬ 
ation processes. As the turbulence inside the cloud decays, 
the cloud loses turb ulent support and collapses inwards (see, 
e.g., iGao et aPl2015h . This is complicated by fragmentation, 
in which dense clumps form and cause the structure of the 
cloud to become highly non-spherical. Since a full model of 
the dynamics of a cloud is beyond the scope of the current 
paper, for the rest of the paper we simply assume u ex t = v eBC - 
Using the power law density field given in Equation [4] 
and the fact that i> 2 sc = 2 GM(< r)/r for M(< r) being the 


mass enclosed inside r, we can compute the escape velocity 
of the cloud as 

2 »iH n, r 2 ~ w 

v esc = 8tt Gn 0 — r 0 -. (13) 

X 6 — W 

An equivalent expression for the ram pressure is 

jj 87 tG / TYln w 1 — w \ ^ /I /(\ 

Pram= 3^(^ nor ° r ) ' (14) 

Applying this to Equation [12] where r = r a t a ii, we can write 


7* stall 


_3_&,\3 f X (3- w) A / Ci \ 
47T q b / V m H SnG J \n 0 rQ J 


(15) 


for a power law density field with gas moving radially in¬ 
wards at the escape velocity, where 9 = 4/(7 — 4ui). 

In Section [7] we compare this model to our simulation 
results to determine whether it holds in a 3D density field. 


3 METHODS 

3.1 Numerical simulations 

We now introduce the numerical simulations designed to test 
the analytic expressions given in the previous section. The 
properties of each of the simulations run are given in Ta¬ 
ble □ In all the simulations in this paper we use a cloud 
of 10 4 solar masses. We leave an analysis of different cloud 
masses to future work, though our analytic results should 
apply to clouds of different masses. All of our simulations 
use a cubic adaptive octree mesh with a coarse spatial reso¬ 
lution 256 3 cells, or 8 levels of the octree, and two additional 
levels of refinement, leading to a total of 10 levels of refine¬ 
ment and a maximum effective spatial resolution of 1024 3 
cells in the most refined cells. We refine based on the Jeans 
criterion, such that if the thermal energy of a cell above the 
maximum level of refinement is lower than its gravitational 
potential, we trigger an additional level of refinement. We 
implement outflow boundary conditions, such that any mat¬ 
ter that leaves the simulation volume is assumed to be lost 
to the system. 


3.2 Initial Conditions 

The setup for our sim ulations is similar to that in 
llffrig fc Hennebehel (l2015ll . In Figure [T] we plot initial hy¬ 
drogen number density profiles for each of the set of initial 
conditions used in this paper. We impose a spherically sym¬ 
metric cloud onto the simulation volume with an inner hy¬ 
drogen number density profile given by a pseudo-isothermal 
sphere n ex t(r,t = 0) = no/(l + (r/r c ) 2 ). Most of the simula¬ 
tions in this paper are run using the “Fiducial” cloud, which 
has no = 9370 atoms/cm 3 and r c = 1.12 pc, with a free-fall 
time of 1.25 Myr, defined as yj 3 ^ p , where p is the mean 
cloud density in the pseudo-isothermal portion of the cloud. 
We impose a cut-off at 3 r c (where n ex t = 0.1 no). The 
temperature inside the inner part of the cloud is set to 10 
K. Outside this region we impose a flat density field at 93.7 
atoms/cm 3 out to 7.6 pc. We then impose a medium outside 
this radius with hydrogen number density 1 atoms/cm 3 at 
a temperature of 4000 K. The total box size is 27 pc, giv¬ 
ing a spatial resolution of 0.1 pc on the coarse grid and a 
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Name 

logi 0 (S»/s 4 ) 

B-field? 

tg/ Myr 

Varying UV Source 

N 00 .B 02 

(no photons) 

/ 

1.25 

N47.B02 

47 

/ 

1.25 

N48.B00 

48 


1.25 

N48.B02 

48 

/ 

1.25 

N49.B02 

49 

/ 

1.25 

Delayed UV Emission 

N48_B02_F2 

48 

/ 

1.25 

N48_B02_F3 

48 

/ 

1.25 


Varying Compactness 


N00_B02_C 

(no photons) 

/ 

0.156 

N48_B02_C 

48 

/ 

0.156 

N00_B02_C2 

(no photons) 

/ 

0.527 

N48_B02_C2 

Infall-Dominated 

48 

/ 

0.527 


N00.B00JF 

(no photons) 

1.25 

N48.B00JF 

48 

1.25 

N49.B00JF 

49 

1.25 

N50.B00JF 

50 

1.25 



Figure 1. Initial hydrogen number density against radius for each 
of the cloud profiles (Fiducial, More Compact and Most Compact) 
used in the simulations. The density field in the Fiducial case ex¬ 
tends to 13.5 pc (the edge of the simulation volume). See Section 
El 2 for details. 


Table 1. Table of the simulations included in this paper. “N” 
refers to the photon emission rate St in all frequency groups in 
photons per second as a power of 10 (with “N00” referring to a 
zero photon emission rate). “B02” refers to a ratio between the 
free-fall time and Alfven crossing time of 0.2, while “BOO” refers to 
a uniform field of 0 //G. “F2” and “F3” refer to simulations where 
the source is turned on at 2 tg and 3 tg respectively (rather than 
after 1 fj as in all other simulations). “IF” refers to an “infall- 
dominated” setup in which there is no turbulence in the initial 
conditions, so the cloud undergoes radial collapse. These runs are 
spherically symmetric due to the absence of turbulence, which 
allows better comparison to our ID analytic models. All simula¬ 
tions are run with the “Fiducial” cloud initial conditions except 
those labelled “C2” (“More Compact”) and “C” (“Most Com¬ 
pact”). The Fiducial, More Compact and Most Compact clouds 
have ratios of sound crossing time t sc t to tg of 0.1, 0.075 and 0.05 
respectively. See Section |3.II for more details. 


maximum resolution of 0.026 pc in the most refined cells. 
The centre of the cloud is located at the centre of the simu¬ 
lation volume. For simulations containing a magnetic field, 
the field is imposed so that the ratio between the free-fall 
time and the Alfven crossing time is 0.2 (denoted in the 
simulation names by B02), corresponding to a value of B = 
20 /.iG at the centre of the sphere and 4.4 /rG at the edge. 
In addition to this, we run one control simulation with no 
magnetic field. A turbulent velocity field is imposed over the 
grid, such that the kinetic energy in turbulence in the cloud 
is approximately equal to the gravitational energy of the 
cloud. The turbulence has a Kolmogorov power spectrum 
with random phases. 

We also run two similar simulations in denser environ¬ 
ments in order to probe the effect of cloud density on the ion¬ 
isation front. The “More Compact” (suffix _C2) and “Most 
Compact” (suffix _C) clouds have ratios of sound crossing 
time to free-fall time of 0.75 and 0.5 of the fiducial setup 
respectivey. This gives free-fall times of 0.75 3 and 0.5 3 of 
the free-fall time in the fiducial cloud (0.527 and 0.156 Myr) 
and cloud radii of 0.75 2 and 0.5 2 respectively, based on the 


expression for the free-fall time given above. The Most Com¬ 
pact cloud has an initial density profile with the same equa¬ 
tion for next as above but with no = 6 x 10 5 atoms/cm 3 
and an edge radius of 3 r c = 0.8 pc. The More Compact 
cloud has no = 5.3 x 10 4 atoms/cm 3 with an edge radius 
of 3 r c = 2.0 pc, with a flat density field of 530 atoms/cm 3 
out to 2.5 pc. As in the fiducial setup we include a diffuse 
medium surrounding the cloud with density 1 atoms/cm 3 
and temperature 4000 K. The More Compact runs have a 
maximum spatial resolution of 0.015 pc and the Most Com¬ 
pact runs 0.0066 pc. 

We run a further four runs using the Fiducial cloud, 
except with no initial velocity field and a magnetic field 
strength set to zero, causing the cloud to be infall-dominated 
and spherically symmetric. These simulations have the suffix 
“IF”. We do this to test our model in a spherically symmet¬ 
ric environment. We include sources of logio(£») = 48,49 
and 50 as well as a run with no source of photons. 

3.3 UV Source Properties 

For each cloud density we run a simulation without photons 
(N00J302, N00_B02_C2 and N00_B02_C). We also run a 
set of simulations with a constant source of photons placed 
at the centre of the simulation volume. For each cloud den¬ 
sity we run a simulation in which an emission rate of ionising 
photons S t = 10 48 s _1 is turned on after one free-fall time 
(N48J302, N48_B02_C2 and N48_B02_C, with t s = 1.25, 
0.527 and 0.156 Myr respectively). This source is modelled 
as a blackbody similar to a BOV star with surface temper- 
ature 33700 K, radiu s 8.3 R© and mass 20 M© (based on 
Isternberg et al J 120031 . calibrated to give roughly 10 48 ion¬ 
ising photons per second). For each photon group we sam¬ 
ple this blackbody and add the appropriate photon emis¬ 
sion rate to each group. In every simulation in this paper 
we assume that the sources have a constant emission rate, 
though massive stars will in practice have a decreasing emis¬ 
sion rate in time as their envelopes expand and decrease in 

















6 S. Geen 


temperature (as we m odel in iGeen et al.1 120151. However, 
ISternberg et al.l ll2003l ) give a more or less constant photon 
emission rate for 3-5 Myr, depending on the lifetime of the 
most massive star in the cluster, which is as long or longer 
than each of the simulations in this paper. 

For the Fiducial cloud we run two further simulations 
identical to N48_B02 but with different photon emission 
rates - a 10 47 photons/s source modelled on a B1V star with 
surface temp erature 285 00 K, radius 4.7 Rg and mass 12 Mg 
(based on lPecaut et al.ll2012f) . and a 10 49 photons/s source 
modelled on 05V star with surface tempe rature 460 00 K, ra¬ 
dius 11.2 Rg and mass 55 Mg (based on lVacca et al.lll996l) . 
These simulations are labelled N47_B02 and N49_B02 re¬ 
spectively. As in simulation N48_B02 we turn on the source 
after one free-fall time in the Fiducial cloud, 1.25 Myr. The 
goal of this is to test the expansion of the ionisation front 
in the cloud for various photon emission rates. In principle 
higher photon emission rates exist for very massive stars 
or more massive clusters, though we find that our 10 49 pho¬ 
tons/s star can easily destroy the cloud and thus any sources 
above this value are expected to behave similarly in runs 
with a turbulent cloud. For the infall-dominated runs, where 
the density of the cloud increases dramatically without tur¬ 
bulent support, we include an extra 10 5 ° ph otons/s source. 
This a blackbody based on a 120 Mg star from lSchaller et ahl 
(119921) with surface temperature 56000 K and radius 22.2 
Rg. 

We also run a simulation identical to N48_B02 but 
without a magnetic field in order to study the effect that 
the held has on the cloud and HII region. We label this 
simulation N48_B00. 

In addition, we run a further two simulations that are 
identical to N48_B02 except that we delay the time the 
source is turned on. These are labelled N48_B02_F2 (where 
the source is turned on at 2 ts = 2.5 Myr) and N48_B02_F3 
(where the source is turned on at 3 ts = 3.75 Myr). We do 
this in order to test whether allowing the cloud to evolve 
further has an impact on the resulting HII region. This is 
motivated by the fact that the cloud structure changes dra¬ 
matically over 3 ts, and we wish to determine whether this 
has an effect on the subsequent expansion of the HII region. 


3.4 Simulation Setup 

The initial conditions are evolved using a modified version 
of the radiative magnetohydrodynam i cs code RAMSES-RT 
llTevssieij 120021 : iFromang et al.ll2006l : iRosdahl et al.l l2013h . 
RAMSES-RT uses a first-order moment method for the ad- 
vection of photons, closing the set of equations with the lo¬ 
cal Ml expression for the radiation pressure tensor. It tracks 
the ionisation states of hydrogen and helium in the gas, and 
couples the interactions between the photons and the gas 
on-the-fly. We split the radiation into three groups, brack¬ 
eted by the ionisation energies of HI, Hel and Hell (13.6, 
24.6 and 54.2 eV for the lower bounds of each), though in 
practice helium ionisation is found to have less of an impact 
than hydrogen ionisation on our results. More energetic pho¬ 
tons do not tend to significantly alter the results since the 
temperature of the ionised gas is determined largely by the 
cooling of photoionised metals in the HII region. In all of the 
simulations in this paper we assume a Solar metallicity ev¬ 
erywhere at all times. We do not consider photons below the 


ionisation energy of hydrogen, nor do we include radiation 
pressure, which we reserve for future work. A reduced speed 
of light of 10 -4 c (= 30 km/s, or 2.4 Ci) is used. We do this 
because the speed of lig ht sets the Courant factor in the 
timestep calculation (see IRosdahl et al.l 120131 . for details), 
and thus a reduced speed of light improves the efficiency of 
our simulations dramatically. We chose the minimum value 
such that the speed of ionisation fronts in our simulations is 
calibrated to be the same as that for a larger speed of light. 

Gas above a hydrogen ionisation fraction of 0.1 is con¬ 
sidered to be pho toionised, and i s set t o a temperature of 
8400 K, following iTremblin et alj (120141 ) for an HII region 
at the same distance from the galactic centre as the Sun. 
This is based on observed HII regions, taking the galactic 
radius of our cloud to be at the same distance from the 
galactic centre as the Sun. For gas below this ionisation 
fraction, non-equilibrium cooling on hydroge n and helium is 
calculated based on the prescription given in IRosdahl et alj 
(120131) . Radiative cooling of metals and background heating 
from the ISM are calculated according to an analytic cooling 
function given in ilffrig fe Hennebella (|2015l ) . This function 
is a combination of the low -temperature cooling function of 
lAudit fe Hennebellel d20 05l) an d the high temperature com- 
ponent of lSutherland fe Dopital (1 1993)4 . giving a cooling and 
heatin g function similar to that used bv ljoung fe Mac Lowl 
l| 20061) . 


EXPANSION IN A SPHERICALLY 
SYMMETRIC COLLAPSING CLOUD 


In the first instance we compare our models to simulations 
in a spherically-symmetric density field without turbulence. 
This is to allow more direct comparison with our ID mod¬ 
els. We do not run the simulations in ID as RAMSES 
does not support spherical coordinates, and thus a direct 
ID comparison between our code and the model is not pos¬ 
sible. RAMSES-RT has taken part in the S TARBENCH 
code comparison project dBisbas et al.l [201~5h and has been 
shown to agree with other codes and analytic theory for 
the expansion of HII regions in a selection of static media, 
in cluding tests that compare our r esults to th e equa tions 


in cluding t ests tnat compare our r esults to th e equa tions 
of Hosokawa fe Inutsukal (120061 4 and iRaga et all (120121 ) . See 
IRosdahl et al.l 1 20131) for additional analytic model compar- 


In order to bridge the gap between these static media 
and a fully turbulent cloud, we compare our equations to a 
source expanding in a cloud without turbulence. It should be 
noted that a cloud with no turbule nt support is c onsidered to 
be unlikely. Observations by, e.g. IPeretto et ahl (120131 ). find 
that clouds tend to be globally virialised. This turbulence- 
free simulation is thus provided largely for the sake of model 
comparison. 

In Simulations N0CLB0CLIF, N48_B00_IF, 
N49_B00_IF and N5CLB0CLIF, we run the Fiducial 
cloud without an initial turbulence field or magnetic field 
but with self gravity (see Section [SJ. In this case the density 
of the cloud incr eases dramat ically on a timescale of the 
freefall time (see iLarsonl Il969l ). We do not provide a de¬ 
tailed algebraic solution for this case. Instead, we compute 
a “non-static” solution to Equation [8] that includes the 
infall velocity of the gas at each radius. We do this by 
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Figure 2. The median radius of the ionisation front over time 
compared to Equation[8]in Simulations N48_B00_IF, N49_B00_IF 
and N50-B00JF. The solid lines show the median radius of the 
ionisation front in each simulation. The dashed lines show the 
“Non-Static Model” solution to Equation [8j which includes an 
external velocity field. 

sampling the time-dependent radial density and velocity 
field in SimulationNOCLBOCLIF and compare that to our 
simulation results. In other words, we use the simulation 
without feedback to provide the density and infall velocity 
for the model. This allows us to test the accuracy of the 
front propagation equation in the absence of a detailed 
theoretical model for the radial collapse of the cloud. 

In Figure [2] we plot this comparison. Our solution to 
Equation [5] (dashed line) follows the evolution of the simu¬ 
lation results (solid line). There is a short time lag between 
our solutions and the simulation results. This is because the 
velocity field of the cloud responds to the HII region, so there 
is some divergence between the simulation without photons 
used to calculate the analytic solution and the simulated 
clouds containing an HII region. In addition, the simulation 
outputs have limited time resolution, and the velocity and 
density fields are linearly interpolated in time. Nonetheless, 
we reproduce the form of the simulation results, in particu¬ 
lar the parabolic arc of the radius over time as the increasing 
density crushes the ionisation front. We also find agreement 
that beyond a certain photon flux the front is able to es¬ 
cape the cloud, as we discuss in more detail in Section [T] for 
the turbulent clouds. In particular, we find that the photon 
emission rate needed to escape the cloud is much higher in 
the infalling cloud than in the turbulent cloud because of 
the dramatic increase in density over time due to the col¬ 
lapse of the cloud, which lacks support against gravity from 
turbulence. 


5 SIMULATION OVERVIEW 

We now review the qualitative properties of simulations in 
the Fiducial cloud and the effect of the magnetic field on the 
structure of the cloud and the HII region. We leave quanti¬ 
tative comparison of the simulation results to our analytic 
models to subsequent sections. 


5.1 Before the First Star 

In this section we briefly review the behaviour of the cloud 
with and without the magnetic field in the absence of pho¬ 
tons as it pertains to the subsequent evolution of the HII 
region. After t = 0, the cloud begins fragmenting under the 
influence of its initial turbulence, magnetic fields and self¬ 
gravity. This fragmentation occurs over ~ ts. In the absence 
of a source of photons, these clumps move towards the centre 
over a further 2 ts- Since the cloud is virialised, this infall is 
due to mass segregation and turbulent decay. The cloud is 
distorted from its initial spherical shape, though maintain¬ 
ing the sharp discontinuity between the cloud gas at above 
100 atoms/cm 3 and the external medium at 1 atoms/cm 3 . 
The shape of the cloud is important insofar as it determines 
how far the HII region has to expand to escape the cloud. At 
tg, we sample the distance from the centre of the cloud to 
the edge along evenly sampled lines of sight (see Appendix 
0 . We find a roughly uniform probability distribution of 
radii from 3 to 12 pc in the Fiducial cloud. Thus the HII 
region must travel at least 3 pc to escape the cloud, and if 
it travels 12 pc it will have completely disrupted the cloud. 

In Figure|3]we plot the cloud structure tg = 1.25 Myr in 
simulations N48_B02 and N48_B00, i.e. the Fiducial cloud 
with and without magnetic fields. At this time no UV pho¬ 
tons have been emitted. In this figure we sample profiles 
along lines of sight as before (Appendix 0 and plot the 
(volume-weighted probability) distribution of density and 
radial velocity at each radius as a colour gradient. The grey 
lines denote the maximum and minimum values at each ra¬ 
dius, and the black line the median value at each radius. The 
red values represent probability bins around the median at 
each 2nd percentile, and the white values the highest and 
lowest percentiles. 

The velocity field is relatively flat in the magnetised 
cloud, whereas in the non-magnetised case flows up to four 
times faster are seen in the non-magnetised cloud. This is 
because, in the absence of magnetic support the cloud must 
rely on support from turbulence. The mean density profiles 
(dashed lines in Figure [3j with and without magnetic fields 
are very similar. Fitting the mean density profiles of both 
clouds at tg to a power law, we find power law indexes (—w 
in Equation [4j of -0.74 in the cloud containing a magnetic 
field and -0.6 in the cloud with no magnetic field. 

However, there are differences in the angular distribu¬ 
tion of matter as shown in the median (solid black line) and 
interquartile ranges (dotted black lines). These drop off more 
quickly in the simulation without magnetic fields. Since the 
mean density profile is the same, it leads us to conclude that 
more mass is found in small clumps in the non-magnetis ed 
clo ud. This agrees wi th the findings of iHennebellel (120131 1 - 
see ISoler et al.l ll2013lf for a discussion. This clumping can 
also be seen in the top panel of Figure [4] though it is clearer 
after the HII region has formed, which we discuss in the next 
section. 


5.2 HII Regions with and without Magnetic 
Fields 

Once the source is turned on at ts, the ionisation front be¬ 
gins to expand. The densest clumps remain embedded while 
the less dense gas is pushed away. There is a competition 
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Figure 3. Probability distribution in density and radial velocity at each radius in simulations with and without a magnetic field. The 
top row shows hydrogen number density, while the bottom row shows the radial velocity, where positive values are flows away from the 
centre. The left-hand figures show the fiducial cloud at % with a magnetic field (simulation N48JB02) while the figures on the right 
are taken at the same time from an identical setup but without a magnetic field (N48-B00), noting that neither simulation contains a 
source of photons before % = 1.25 Myr. In each plot we sample the density and radial velocity along evenly spaced lines of sight from 
the source position (see Appendix [Dt. The grey lines denote the maximum and minimum values for all lines of sight at each radius. The 
solid black line shows the (volume-weighted) median value at each radius. The dashed black line is the mean value at each radius. The 
dotted black lines show the values at the 25th and 75th percentile. The shaded region goes from white (maximum or minimum at each 
radius) to red (median). 


bet ween the acceleration of the clumps by the rocket effect 
(see lOort fe Spitzeiiri955l '). in which UV photoevaporation 
from the surface closest to the source causes the clump to 
accelerate away, and the effects of gravity. The most mas¬ 
sive clumps, unless resisted by the UV photons, will tend 
to m ove further inwards due to mass segregation (ISpitzed 

1 1964 . 

The increased fragmentation of the cloud when no mag¬ 
netic field is present provides more channels of low-density 
gas between the dense clumps. This in turn allows the HII 
region to escape preferentially through these channels. This 
can be seen in the bottom panels of Figure 2] where we plot 
the maximum extent of the HII region along the line of sight 
of the projection image as a cyan (light blue) contour. In the 
simulation with a magnetic field, by contrast, the smoother 
density field causes the HII region to become (relatively) 
more spherical. 


For most of the simulations in this paper we include a 
magnetic field. While there are major qualitative differences 
between the results with and without magnetic fields, the 
mean density field in each simulation is similar. As a result 
quantities such as median ionisation front radius, mass of 
ionised gas and momentum of the cloud gas are broadly 
similar. We discuss comparisons between these quantities 
and analytic theory in the next section. 


6 INFLUENCE OF PHOTON EMISSION RATE 
AND CLOUD COMPACTNESS 

In this section we compare our simulations to our analytic 
models, focussing on the influence of varying photon emis¬ 
sion rates and cloud compactness on the properties of HII 
regions. In the first instance we compare our simulations to 
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Figure 4. Sequence of projections of the “fiducial” cloud showing maximum hydrogen number density along the line of sight. The cyan 
(light blue) contour marks the edge of the HII region (measured as a hydrogen ionisation fraction above 0.1) when projected onto the 
image. A red circle indicates the position of the UV source if one is included in the simulation. The left column is simulation N00_B02, 
the middle column is N48J302, and the right-hand column is N48-B00. The top row shows each simulation at 1.25 Myr (1 fg) and the 
bottom column at 3.25 (tg + 2 Myr). All images show the full cubic simulation volume of length 27 pc. The presence of a source of UV 
photons dramatically alters the overall structure of the cloud, while the magnetic field alters the filamentary structure and the shape of 
the HII region. 


analytic models that assume a spherically symmetric power 
law density field with no velocity or pressure terms outside 
the ionisation front. The point of this exercise is to deter¬ 
mine to what extent previous models for the expansion of 
ionisation fronts hold in more complex media. In subsequent 
sections we discuss at what point these models break down 
and introduce new ones that fit the simulations better. We 
invoke Equations [5] to Q given in Section \2 . 2 1 

We fit our simulated clouds to a power law density field 
(Equation [4]) using the spherically-averaged density profile 
sampled from the simulation output at the time the source is 
turned on, sampling inside a sphere of radius 13.5 pc around 
the source (50% of the distance to the edge of the box). 
These fits are not perfect descriptions of the actual density 
field, but they give a reasonable agreement to the density 
field while allowing comparison to analytic models. It also 
allo ws a p oint of c ompa rison with observed clouds, though 
as lTremblin et akl (120141 ) note, there exists some degeneracy 
with regards to modelling the age and radius of observed 
HII regions. We do not directly address comparsion with 
observations in this work. In the fiducual cloud at tg we 
find a power law index of -0.74. This gives a solution to 
Equation [5] where n oc S°' 18 t 0 ' 72 . 


Our simulations include a static source of UV pho¬ 
tons in the centre of the simulation volume, turned on after 
one cloud free-fall time. In reality, the UV emission rate 
will be determined by the mass of the most massive star 
dVacca et al.l 1 19961 : iMartins et al.l 120051) , which in turn de¬ 
pends on the mass of the cluster formed as well as how the 
IMF is sampled. We will treat self-consistent star forma¬ 
tion in future work. However, it should be noted that as 
the most massive stars fall towards the centre of the cloud 
as the cluster undergoes mass segregation, it is a reason¬ 
able approximation to assume that the majority of the UV 
emission rate is coming from the centre of the cloud. 

In Figure [5] we show the density field of the simu¬ 
lation for different timesteps in the simulations N00_B02, 
N47_B02, N48_B02, and N49_B02 - that is, for simulations 
in the Fiducial cloud but with varying source emission rates 
S* = 0, 10 47 , 10 48 , and 10 49 s _1 . The first and third of 
these simulations are repeated from Figure [I] For the weak¬ 
est source (10 47 s _1 ), the ionisation front is insufficient to 
resist the infalling clumps in the right hand side of the fig¬ 
ure. As a result the HII region expands in one direction only 
as a “blister” region, while for larger UV emission rates, the 
ionisation front expands in all directions. The most diffuse 
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gas is pushed away as a dense shell, while the most massive 
clumps remain embedded inside the HII region. After 3 Myr 
in the 10 49 s _1 simulation, the cloud is nearly entirely de¬ 
stroyed save for a few cometary clouds, while the 10 47 s _1 
source has barely changed the structure of the cloud when 
compared to the run with no UV source. 


6.1 Modelling the Ionisation Front Radius 

We plot the properties of the HII region with varying pho¬ 
ton emission rate in Figure [6] In Figure [5] we have already 
seen that there is a considerable scatter in the radius of the 
ionisation front with angle. As a result there is no single 
ionisation front radius as in the ID models, but a distri¬ 
bution of radii. Therefore, we plot the median radius of the 
ionisation front, measured by sampling the radius of the ion¬ 
isation front along several lines of sight. We use the median 
rather than the mean as the latter biases towards extreme 
values, causing our results to overestimate the radius of the 
ionisation front. Rather, we are interested in whether the 
front is able to escape over the majority of lines of sight. 
Comparing these results to the solution of Equation [5] given 
the power law fit to each simulation as described above, 
we find a reasonable agreement between the analytic theory 
and the simulation results at early times in the runs with 
10 48 and 10 49 photons/s sources. A key discrepancy between 
the power law density fit and the 3D cloud structure is the 
presence of a sharp discontinuity the at cloud edge between 
the dense cloud gas and the diffuse external medium. This 
discontinuity lies at a radius varying from 3 to 12 pc de¬ 
pending on the line of sight from the source. As a result the 
ionisation front expands rapidly into the external medium 
once it reaches this edge, which happens at around 1 Myr 
with the 10 49 s _1 source and 1.5 Myr for the 10 48 s _1 source. 
We discuss a model that corrects for this in Appendix m 
though for reasons of simplicity we do not use this model in 
our analysis. 

The results for the simulation without a magnetic 
field (N48_B00) are very similar to the results with one 
(N48J302) during the early expansion phase, since the mean 
density profiles are similar. However, as we discussed in Sec¬ 
tion [5T] the density held is more fragmented without a mag¬ 
netic field. This means that the ionisation front can leave the 
cloud more rapidly. As a result, it accelerates away from the 
power law solution sooner than in the run with a magnetic 
held. 

The radial evolution of the 10 47 photons/s run in the 
top panel of Figure[I>]is significantly different from the power 
law model. While the ionisation front leaves the box over 
a small fraction of the solid angle around the source, the 
median radius of the ionisation front stagnates. In other 
words, for at least half of the solid angle around the source, 
the ionisation front remains trapped by the cloud, at least 
for 4 Myr after the source is turned on. We return to this 
issue in Section 0 

6.2 Ionised Mass and Momentum 

We can also compare the ionised mass and momentum added 
to the system by the HII region to the analytic expressions 
given in Equations [6] and [7] It is worth repeating that even 


for an ionisation front that expands to engulf the entire 
cloud, most of the mass in the cloud is found the neutral 
shell. This is because, from Equation [3] m = no(r B /n) 3 , 
and hence most of the mass displaced by the ionisation front 
is pushed onto the dense shell around the cloud. As with 
the radius of the ionisation front, the mass in ionised gas for 
S'* = 10 47 s _1 is overestimated, since the ionisation front is 
trapped in the direction where most of the mass in the cloud 
is found. The non-zero mass in ionised gas in N48_B00 at 
t=0 is because the gas in the diffuse phase is hot enough to 
become collisionally ionised. This does not affect our subse¬ 
quent results for the photoionisation of the cloud. 

For the momentum in the cloud we find a good fit be¬ 
tween the analytic model and the simulation results un- 
til mass begins to leave the simulation volume. As in 
llffrig fe Hennebehel d2015h we give momentum as the total 
momentum in radial flows in the simulation. The momen¬ 
tum in the 10 49 s _1 simulation approaches the momentum 
found when a supernova goes off in the cloud (see again 
llffrig fe Hennebellell2015h . The 10 47 s _1 source is insufficient 
to displace a significant quantity of mass and the momen¬ 
tum of the cloud is not very different from the case without 
a source of photons at all. The results of the simulations 
with sources of 0 and 10 4 ' photons/s are not identical due 
to the nonlinear nature of the flows in the cloud, but they 
both fluctuate between 2 and 4x 10 42 g cm/s of momentum. 
The momentum in simulations N48_B00 and N48_B02 is 
very similar until the ionisation front leaves the simulation 
volume, which it does earlier in run N48_B00 as explained 
in the previous section. 

In Figure [7] we plot the momentum and ionised mass 
in simulations with increasingly compact initial conditions. 
Varying the cloud compactness gives a similar effect to vary¬ 
ing the photon emission rate. The mass in ionised gas is well 
captured by the Power Law model for the Fiducial and More 
Compact clouds, whereas in the Most Compact case we find 
a negligible quantity of ionised gas. Similarly, for the mo¬ 
mentum, only in the Fiducial case does the Power Law model 
accurately predict the amount of momentum added to the 
system by the HII region. In the More Compact and Most 
Compact clouds, the momentum in the system is dominated 
by the fluctuating momentum of the cloud itself. In simula¬ 
tion with 10 48 photons/s in the More Compact cloud, there 
is an initial divergence in the results but this only maintains 
the momentum in radial flows, compared to a predicted drop 
in momentum in radial flows when there is no source of UV 
photons. It should be noted that the Power Law model does 
not predict a large increase in momentum in the Most Com¬ 
pact cloud, and that the model lies within the variation in 
momentum over time in this cloud. We return to reasons 
why the HII region is ineffective at driving flows in compact 
clouds in the following section. 

6.3 Delayed UV Emission 

One additional test we perform is to determine what hap¬ 
pens when we allow the cloud to relax over a longer period 
of time before we turn the source on. To do this we use the 
Fiducial cloud with a source of S* = 10 48 s _1 . In addition to 
N49_B02 we run a further two simulations in which we turn 
on the source at 2 and 3 tg (2.5 and 3.75 Myr), N48_B02_F2 
and N48_B02_F3 respectively. At 2 tg the power law index 
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of the cloud density profile becomes -0.87, and at 3 tg it 
becomes -1.88 (compared to -0.74 at tg). Nonetheless, the 
resulting expansion of the HII region does not significantly 
change between these simulations. There is a short period 
during which the ionisation front stalls in both delayed sim¬ 
ulations, but after 1 Myr the front breaks out and expands 
similarly to N48_B02. This suggests that the properties of 
a HII region are not very sensitive to when stars are formed 
in a cloud with a given initial structure, though there may 
be a short period during which the front must escape the 
increased density in the centre of the cloud. 


7 STALLED EXPANSION 

In the previous section we compared our models assuming 
a static power law density profile to our simulations in the 
Fiducial cloud. These models provided a good fit to the sim¬ 
ulations provided the UV photon emission rate was above 
10 48 photons/s. However, the front around the 10 4 ' pho¬ 
tons/s source stalled and was unable to expand further for 
4 Myr. In this section we compare our model for expansion 
in a medium with a non-zero external velocity field given in 
Section m in order to explain this effect. 

In Figures [ 6 ] and [ 7 ] we plotted the Power Law model as¬ 
suming no external pressure terms. This model fit the sim¬ 
ulation results well for large photon emission rates or less 
compact clouds, but not so well for weak sources or dense 
clouds. In this section we address this by solving Equation 
[8] including a velocity field. We plot this solution as a “Non- 
Static Model” (as opposed to the model in the previous sec¬ 
tion that assumed a static cloud with no velocity field) next 
to our simulation results in Figure [8] 

In each solution to this equation use two approxima¬ 
tions for the density and velocity field. In the “Power Law” 
model, we assume a power law density profile as in the previ¬ 
ous section, with a velocity profile calculated as the escape 
velocity at each radius as in Equation [13] This model is 
highly idealised, and so we also plot the “Sampled” model, 
which samples the time-dependent density and velocity field 
from an identical simulation without a source of photons. 
The values sampled are taken from the spherically-averaged 
radial density and velocity profiles in each output and inter¬ 
polated in radius and time. The “Power Law” model pro¬ 
vides a close comparison with our simple analytic model, 
whereas the “Sampled” model gives a closer match with our 
simulated cloud. Both models however assume a spherical 
density and velocity field. We discuss the limitations of this 
assumption below. 


7.1 Free vs Stalled Expansion 

In the top panels of Figure[8]we use the Fiducial cloud while 
varying the photon emission rate. The 10 48 and 10 49 pho¬ 
tons/s sources cause the ionisation front to break out of the 
cloud, while the front stalls with a 10 47 photon/s source. 
Similarly, in the bottom panels, the More Compact and 
Fiducial cloud allow the ionisation front to expand freely 
while in the Most Compact cloud the front stalls. Our so¬ 
lution to Equation [8] including a velocity field thus matches 
our simulation results better for the weak source and most 
compact cloud. 
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Figure 6. Properties of the HII region over time for varying UV 
photon emission rates, with the same simulations as in Figure [5] 
plus N48_B00. In the top panel is the median radius of the ioni¬ 
sation front across randomly sampled lines of sight. In the middle 
panel is the mass in ionised gas. On the bottom is the momentum 
in radial flows in the simulation. Each colour corresponds to a dif¬ 
ferent photon emission rate. The simulation results are shown as 
a solid line. We plot as a dashed line the Power Law model (see 
Section l2.2> for each photon emission rate. The time given is from 
the time the source is turned on (tg = 1.25 Myr after the start 
of the simulation). We overplot as dotted lines the time in each 
simulation where the ionisation front first leaves the simulation 
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Figure 5. As in Figure [T] but for varying UV photon fluxes. From left to right the columns are for simulations N0CLB02, N47-B02, 
N48-B02 and N49_B02. The rows show each simulation at tg + [0,1,2,3] Myr respectively. The cloud at iff is shown in the top left panel 
in Figure [4] 


The Power Law and Sampled models both predict the 
radius at which the ionisation front stalls in the simula¬ 
tions where this occurs. The Sampled model also reproduces 
the expansion and contraction of the ionisation front in the 
case where it stalls. We discuss the cause of this in Section 
El The Sampled model predicts incorrectly that N48_B02 
should stall for around 2 Myr, since it assumes that that the 
density field is spherically symmetric, meaning that the ion¬ 
isation front cannot escape through channels of low density 
as it does in the simulation. 

The Power Law model tends towards the value for 
v stall given in Equation [15] In this equation, r sta ii varies 
as (5y 4 /no)^. t/> i s larger for steeper density profiles. This 
means that the effect of cloud density and photon emission 
rate is slightly enhanced in the Most Compact cloud, which 
has w ~ 1 (ip = 0.8), compared to the Fiducial cloud, which 
has w ~ 3/4 ('!/’ = 0.72). We give the value for r s t a ii in each 
simulation in Table [2] In none of our simulations do we find 
that the escape velocity exceeds the sound speed of the gas, 
and hence in all simulations r s t a ii > r s . 

In the following subsections we discuss the reason why 
the ionisation fronts in some simulations seem to expand 


Simulation 

^cloud/ P c 

r sta ll/pc 

^stall/^cloud 

N47.B02 

3.0 

3.55 

1.18 

N48.B02 

3.0 

6.24 

2.08 

N49.B02 

3.0 

10.4 

3.47 

N48.B02.C2 

1.7 

2.83 

1.68 

N48.B02.C 

1.0 

0.433 

0.577 


Table 2. Table comparing the minimum cloud radius r c f ou< f (as 
measured at tg in each simulation) to the radius at which the 
ionisation front stalls in the analytic model r s t a ii. If the cloud 
radius is smaller than the stalling radius, the front can escape the 
cloud. 

freely into the external medium, while in others it stalls. We 
also discuss why r s t a ii in the simulations and in the Sam¬ 
pled model appears to oscillate (causing the HII region to 
“flicker”). 

7.2 Role of the Diffuse Medium 

The assumption of spherical symmetry in the density and 
velocity field breaks down when we reach the cloud edge at 







Photoionisation Feedback in Molecular Clouds 13 




Time / Myr Time / Myr 




Figure 8. The median radius of the ionisation front over time compared to Equation [8] The top panels show simulations N47J302, 
N48-B02 and N49-B02. The bottom panels show simulations N48-B02, N48_B02_C and N48_B02_C2. The left panel uses a model 
assuming a power law density profile and velocity profile given by the escape velocity at each radius r)/r. The right panel 

models the density and velocity field of each simulations directly by sampling the density and velocity at each radius and time in identical 
simulation without a source of photons. The solid lines show the median radius of the ionisation front in each simulation. The dashed 
lines show the “Non-Static Model”, in which we solve Equation [8] including an external velocity field. 


v cloud ■ The cloud edge is a sharp transition of a fraction of a 
parsec between the dense cloud gas and the diffuse external 
medium. The Power Law model does not resolve this edge 
because r c ioud varies over a large range (3 to 12 pc in the 
Fiducial cloud). If r sta ii/r c i ou d > 1, the front will enter the 
diffuse medium, where the ram pressure becomes very low 
since the density and velocity are much lower than inside 
the cloud. Hence beyond this radius the front will no longer 
stall, but expand very rapidly outwards. 

In Table [2] we give the values for r c ioud, fstaii and 
v stall / r cloud • The latter value is close to 1 in simulation 
N47_B02. This means that (on average) the front is only 
just expected to escape the cloud. If we compare ?' a t a ii in this 
simulation to Figure [SI we find that it has not yet reached 
v stall■ However, in Figure^ we do see a sudden expansion of 
the mean ionisation front radius at 4 Myr, which can be at¬ 
tributed to the ionisation front breaking out of the cloud. In 
simulation N48_B02_C, the ratio is 0.577. Thus the ionisa¬ 
tion front is expected to remain trapped inside the cloud. In 
all other simulations r s taii is much larger than r c ioud. This 


means that the stalling criterion does not affect the ioni¬ 
sation front beyond this radius and expand freely into the 
diffuse medium. 

We thus have a criterion to determine whether the me¬ 
dian radius of the ionisation front should follow the static 
Power Law model or a model that includes ram pressure 
from the external velocity field. This criterion is very impor¬ 
tant in determining the other properties of the HII region 
and cloud, such as momentum injected into the ISM as well 
as the strength of star formation feedback. In the next sec¬ 
tion we discuss the role that density inhomogeneities play 
in regulating the expansion of the HII region. 

7.3 Role of Dense Clumps 

Our simulations (and observed molecular clouds) are highly 
non-spherical, with filamentary and clumpy structures. 
While our spherically symmetric ID analytic models work 
well despite this, there are some effects caused by these in¬ 
homogeneities that we now discuss. 
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Time / Myr 

Figure 7. Properties of the HII region over time for varying 
cloud compactness, as in Figure [6] Each colour corresponds to 
a different cloud “compactness”. The solid lines show the simu¬ 
lation results for each cloud compactness with a photon source. 
The dotted lines show the simulation results for the same cloud 
without a photon source. The dashed lines show the Power Law 
model (see Section l2.2l) for each cloud compactness. The time axis 
begins at tg after the start of the simulation, for tg as defined in 
Table |H 

Dense clumps can resist the expansion of the HII re¬ 
gion due to having a higher density than the rest of the 
cloud and a smaller surface area. The most massive clumps 
remain embedded inside the HII region as cometary clouds. 
They block the expansion of the HII region over the solid 
angle subtended by the clump from the source, creating a 
“shadow” of neutral gas behind them. If this angle is small 
enough, it does not affect the median radius. However, if 
they move close to the source they can increase the density 
into which the front expands (n ex t), especially if they pass 
through the position of the source. This prevents the ioni¬ 
sation front from ionising and expelling a larger fraction of 
the cloud. 

In Figure [9] we plot the time evolution of the distance 
between the source and the edge of the nearest clump. We 
define this as the smallest distance between the source and 
a cell above a given density threshold. This threshold is 
10 6 atoms/cm 3 in the Fiducial cloud, 10 7 atoms/cm 3 in 


the More Compact cloud, and 10 s atoms/cm 3 in the most 
compact cloud. This is because the cell sizes in the denser 
simulations are smaller, allowing us to resolve gas at higher 
densities. The thresholds for defining which gas is in dense 
clumps are found such that there is a 1% chance that a ray 
cast from the source to the edge of the cloud at tg in a ran¬ 
dom direction will encounter a cell of at least this density. 
This is seen in the right hand plot of the same figure, where 
we plot the probability distribution function of maximum 
density along the lines of sight from the source position to 
the edge of the simulation volume. The reason the ionisation 
front is non-spherical is due to this probability distribution 
function of densities around the source. 

As we increase the photon emission rate, the clumps are 
accelerated further from the source, meaning that they play 
little role in the expansion of the ionisation front in other 
directions. As such, for higher photon fluxes in the Fidu¬ 
cial cloud, the front can expand freely in most directions. 
The 10 4 ' photons/s source is insufficient to push away the 
clumps. The radius at which the front stalls is thus affected 
in part by the motion of the clumps. The radius of the ion¬ 
isation front oscillates in Figure [8] as the clumps orbit close 
to the source. 

In the denser clouds, the source is largely unsuccessful 
at significantly altering the trajectory of the dense clumps. 
In the More Compact cloud, the clumps remain far enough 
from the source for the ionisation front to escape over at 
least half of the solid angle around the source. By contrast, 
in the Most Compact cloud the clumps orbit close to the 
source, passing through it multiple times. As a result the 
source flickers on and off, as seen in the simulations results 
and the “Non-Static” model in Figure [8] 

The motion of dense clumps inside the cloud thus plays 
a role in determining the shape and behaviour of the HII 
region. However, the broad behaviour of the HII region can 
be described well with a spherically symmetric model that 
takes into account the distance between the source and the 
edge of the cloud. 


8 EVAPORATION OF STAR-FORMING 
CLUMPS 

Molecular clouds are the sites of star formation in galaxies. 
An important role of UV photoionisation is thus the self¬ 
regulation of star formation. This introduces feedback cycles 
from stars, which both photoionise the clouds in which they 
formed and, if the ionisation fronts escape the cloud, affect 
the evolution of nearby clouds. This feedback can either be 
positive, as in the case of shock compression enhancing the 
star formation rate, or negative, as in the case of UV photons 
evaporating clouds that would otherwise be star forming. In 
this work we do not simulate star formation directly, nor do 
we place sources of UV photons onto star particles. We thus 
do not reproduce these feedback cycles directly in our sim¬ 
ulations. Instead we discuss the potential for star formation 
in our simulations based on stability criteria, and compare 
this to analytic models from the literature. 
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Figure 9. Left: Minimum distance from the edge of the nearest dense gas clump to the source position in simulations over time. Right: 
the probability that a ray cast from the source in a random direction at iff in each simulation will encounter material of at least a given 
density in each of the clouds of a given compactness. 


8.1 Jeans Unstable Mass 

Our highest resolution in the Fiducial cloud is 0.026 pc (= 
5 500 AU), which is bel ow th e radius of star-forming cores 
(IWard-Thompson et al.fl994|j but insufficient to follow their 
collapse into protostellar cores. We hence argue that any gas 
cell that is Jeans unstable is to be considered potentially 
“star-forming”. The Jeans stability criterion states that if 
the free-fall time of any given part of the cloud is smaller 
than the sound-crossing time, the cloud is vulnerable to frag¬ 
mentation and collapse. We set the Jeans Length Aj to our 
cell length and count the mass in cells for which 

< 1 ( 16 ) 

where c s is the sound speed in the cell and p is the density. 
For example, a gas cell at our maximum resolution in the 
Fiducial cloud with temperature 50 K and 2x10 s atoms/cm 3 
is at the Jeans stability limit. In the More Compact and 
Most Compact clouds the size of a cell at the maximum spa¬ 
tial resolution is reduced by a factor 0.75 2 and 0.5 2 respec¬ 
tively, proportionally to the radius of the cloud. It is impor¬ 
tant to note that not all mass that fulfills the J e ans criterion 
will end up in stars. iMatzner fe McKee] (120001 ): lAlves et all 
(2007) suggest that the star-formation efficiency (SFE) of 
cores is around 30%. The mass in Jeans-unstable gas is thus 
an overestimate by a factor of 3-4. 

We plot the total mass in Jeans unstable gas in Fig¬ 
ure EH In the Fiducial cloud with varying photon emis¬ 
sion rate, the amount of potentially star-forming mass drops 
with increasing photon emission rate, even with the 10 4 ' 
s ' 1 source, whose median ionisation front radius is trapped 
by the cloud. Mass in dense clumps is evaporated and re¬ 
distributed to the shell around the ionisation front. We do 
not find that this shell is unstable to fragmentation (see Ap- 
pendix[C|. In cases where the ionisation front leaves the sim¬ 
ulation volume, a certain amount of mass disappears from 
the simulation, as in the 10 49 s _1 source. However, many of 
the most massive clumps remain embedded inside the HII 
region (see Figure [5]). 


Aj \JGp 


In the More Compact and Most Compact clouds, the 
radiation has little effect on the Jeans unstable mass. Even in 
the case of the More Compact cloud, where the front escapes 
the cloud, the densest gas in the simulation remains largely 
unaffected by the UV photons. In the previous section we 
explained how, despite the fact that the ionisation front can 
escape the More Compact cloud, the trajectories of dense 
clumps in the cloud are largely unaffected by the HII region. 


8.2 Modelling Clump Evaporation 


In this subsection w e com p are our simulation result s 
to iBerto ldi fe McKeel (fl990h : iLefloch fe Lazarefj (1994); 
Ijohnstone et al. ( 1 998lf : fwh i t. worth fc Zinnecken ( 2004 ), 
who give analytic models for the evaporation of dense gas 
clouds by sources of UV radiation. In our simulations we 
do not track individual clumps of gas. This requires special 
treatment in Eulerian codes such as RAMSES, which do 
not trace fluid parcels as Lagrangian codes do. Instead, we 
sample the dense clumps in our simulation in one snapshot 
at is and model the i r mas s evolution using the equations 
of lBertoldi fe McKee! (Il990l f . We define our dense clumps to 
be t hose portions of gas in th e cloud that are Jeans unsta¬ 
ble. [Bertoidr^McKeg (ll990Tl give, for thermally supported 
clumps, a clumps mass with a fraction 


(1 - f/fev) 5/3 


(17) 


of its mass at t = 0 . t ev , the characteristic evaporation time 
of the clump, is given by 0.448#iCg R^ 5 m^ 5 Myr, 

where 9 1 is a factor close to unity, cs is the sound speed of 
the clump in km/s, S 49 = S*/10 49 , Ri is the distance of the 
clump to the source in pc and mi is the initial mass of the 
clump in Mq. For the clumps in our simulation, i ev varies 
from between 1 to 30 Myr. 

We sample the clumps in our simulations at tg using 
the following clump finding algorithm. We first select all gas 
cells above the Jeans stability threshold as given previously. 
We then assign all contiguous cells above this threshold to a 
single clump. For each clump we calculate the total mass of 






































16 S. Geen 


the cells inside it, its radius (measured as the maximum dis¬ 
tance from the densest cell in the clump to any other cell) 
and the distance between the densest cell and the photon 
source. We find that at tg the clump mass mi, radius n 
and distance to source Ri are related as mi oc r*' /2 (i.e. a 
Larson relation) and mi oc R 1 (measured from our simu¬ 
lation), though there is a large degree of scatter in the latter 
relatio nship. This is in agreement with iBertoldi fc McKeel 
l|l990h . 

From tff, we allow the mass of each clump to vary ac¬ 
cording to the following model. First, we assume that the 
mass of the clumps increase due to accretion as m(t) = 
mi(l + t/tacc), where face is set to 4.7 Myr (the free-fall time 
for a density of 93 atoms/cm 3 , a sphere of which surrounds 
our cloud in the initial conditions). We omit this accretion in 
the Most Compact cloud since there is almost no gas in this 
phase in our initial conditions. Secondly, as the ionisation 
front reaches the clump (using the Power Law model for the 
expansion of the ionisation front) we turn off the accretion 
and begin evaporating the clump according to Equation 1 171 
We do not attempt to follow the orbits of the clumps as 
this requires more detailed modelling of the cloud dynamics 
that is beyond the scope of this paper. Instead, we make 
the simplifying assumption that the clumps remain at the 
same distance from the source at all times. We then sum the 
masses of each clump at each time and overplot as a dotted 
line in Figure flOl 

For the Fiducial cloud in the left-hand panel of Figure 
EB this model is reasonably successful at reproducing the 
mass of Jeans Unstable gas measured in our simulations. 
It overestimates the mass loss for the 10 4 ' s -1 since the 
front stalls and is unable to ionise clumps in the direction 
of the densest parts of the cloud. We overpredict the Jeans 
Unstable mass in the 10 49 s _1 simulation after a few Myr. 
This is because the HII region expels some massive clumps 
from the simulation volume, meaning we are no longer able 
to track this mass. The more clumpy nature of the cloud 
without a magnetic field causes it to produce more Jeans 
unstable mass than the simulation with a magnetic field, as 
the cloud fragments more, while the ionisation front escapes 
preferentially in directions with lower densities. 

For the More Compact and Most Compact clouds, our 
models overpredict the amount of mass evaporated. For the 
Most Compact cloud, the model with no photons (the upper 
light blue dotted line) fits better than the model with pho¬ 
tons (the lower light blue dotted line). In the More Compact 
cloud, both models fail to reproduce the simulation results. 
As we find in the previous section, while the ionisation front 
in the More Compact cloud escapes the cloud, it does not 
affect the position of the dense clumps in the cloud. Thus 
our spherically-averaged model for the ionisation front ra¬ 
dius overestimates the effect that the ionisation front has on 
the clumps inside the cloud. 

Our measurements of the Jeans unstable mass in our 
simulations predict that the main effect of HII regions 
in clouds should be to evaporate star-forming clumps, re¬ 
ducing the star formation efficiency of the cloud. These 
measurements are backe d up by a simple model based on 
IBertoldi fe McKej dl990h . However, in cases where the ioni¬ 
sation front stalls this simple model (which does not include 
the stalling radius) overpredicts the evaporation of dense 
clumps in the cloud. 


9 DISCUSSION 

We now discuss the further consequences of our results and 
some limitations to our work. While our setup is deliber¬ 
ately simplified in various ways in order to make analytic 
comparisons possible, there are various aspects of feedback 
in molecular clouds that we omit that nonetheless are ex¬ 
pected to play an important role. 

One important effect that we leave out of our simulation 
is self-consistent star formation. Instead, we model emission 
from stars by assuming that a single black body source of 
photons is turned on in the centre of the cloud after one free- 
fall time. The most massive stars in a young cluster tend to 
be found in_the cent re of the cluster thanks to mass segre¬ 
gation (|Spitzerlll969 'l. and so it is not entirely inconsistent 
to use the approach given in this paper provided that the 
time that the most massive objects fall into the centre of the 
cloud is shorter than the time the ionisation front takes to 
escape the cloud. In practice, this criterion is best fit by the 
densest cloud, where the crossing time is shortest. 

However, in placing our source of photons by hand we 
miss three important points. The first is that we miss the 
early phase of star formation inside the clumps, in which 
the properties o f the stars are sensitive to the radiative flux 
fe.g. lKetdl2003) . 

The second is that the emission rate of UV photons is 
roughly proportional to th e cube of the mas s of the most 
massive star in the cluster (lYacca et alJll99~3 ). and thus we 
are unable to comment on the link between the cloud prop¬ 
erties and the amount of feedback from UV photons in the 
cloud’s star formation cycle. Observed ultracom p act HII re¬ 
gions are expected to b e short-lived llWood fe Churchwelll 
1 19891 : I Walsh et all [l995). Assuming observed clouds can 
reach the densities found in our compact cloud, one pos¬ 
sibility is that the UV emission rate of the cluster formed in 
the cloud increases due to star formation until it is sufficient 
to break out of the cloud. This scenario requires simulations 
with self-consistent star formation to test, which we leave 
for future work. 

The third is that our source is not gravitationally bound 
to the cloud, and this means that we overestimate the extent 
to which the stars are separated from the dense clumps in 
which they form. In addition to this, our cloud is modelled 
as an isolated sphere, and so star-forming clumps accrete 
only o n the order of the f ree-fal l time in the dense gas in the 
cloud. iDobbs fe Pringlel (l2013ll suggest that clouds should 
accrete as they pass through the spiral arms of galaxies. 

In order to compensate for the lack of self-consistent 
star formation in our simulations, we locate regions of our 
simulations where there is Jeans unstable gas to predict 
where stars should form, and how much potentially star¬ 
forming mass should be lost by evaporation. We compare 
this to a simpl e model for the evaporat ion of dense gas 
clumps based on IBertoldi fc McKeel (119901 1. When analysing 
the response of the Jeans unstable mass in our simulation 
to UV photoionisation, we do not find examples where the 
mass of unstable gas is increased by the presence of a UV 
source, and the feedback from UV photoionisation appears 
to be entirely negative (i.e. it reduces the mass available to 
form stars). We cannot discount a temporary boost in star 
formation r ates from the initia l compre ssion of clouds a s 
predicted bv lWalch et al.l ll2012ll . though, bale et all d2015l l 
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Figure 10. The total mass of Jeans-unstable gas over time for each simulation, starting from the initial conditions. The left-hand 
panel shows simulations in the case where the photon emission rate is changed, N00J302, N47J302, N48J300, N48J302 and N49J302. 
The right-hand panel shows simulations in which the cloud compactness is changed, N0CLB02, N48-B02, N48_B02_C2 and N48_B02_C. 
Overplotted as dotted lines is the model for each UV photon emission rate as described in Section [8] 


argues that we should be cautious in how we define “trig¬ 
gered” star formation. Similarly, we do not find that the 
compression of the cloud into a shell creates Jeans unstable 
gas, though this i s still a subject of ongoing debate (see, e.g. 
lKin^&I<hn|201^)jElmegreen (ll994l l; lHosokawa fe Inutsukal 
( 20061') : Tremblin et al.l 1 2014 1 predict that under certain 
regimes the shells around ionisation fronts should become 
unstable to fragmentation, though we do not find that this 
does not occur in our simulations. Even if the shell does 
not become unstable, the momentum transferred to the ISM 
from the expansion of the ionisation front can be anything 
up to that injected by a supernova, helping to maintain tur¬ 
bulence in the ISM I Gritschneder et alJl2009h . This in turn 
could help regulate star formation on a global ISM scale 
over timescales longer than the life of the cloud, even if the 
influence of stellar feedback from UV photoionsation inside 
a cloud has a net negative effect. 

The results of our model for the stalled expansion of 
ionisation fronts have some important consequences for star 
formation. Firstly, if ionisation fronts are unable to dis¬ 
rupt the cloud in which they form, they will be unable 
to prevent further star formation, at least until the point 
at which the emission rate from UV photons is sufficiently 
large to overcome the infall velocity. This is one possible 
explanation for t h e presenc e of observed “super star clus¬ 
ters” (e.g. IKeto et al.1 1200511 - if massive, compact clouds 
are formed in starburst galaxies, photoionisation feedback 
will be unable to significantly disperse these clouds before 
they form the bulk of their stars. Secondly, the environment 
into which massive stars explode as supernovae changes the 
properties of the resulting blastwave. In circumstellar envi¬ 
ronments pre-processed by UV photoio nisation, more en ¬ 
ergy is retained by the blastwave (e.g. Rogers & Pittarc| 
120131: iGeen et al.ll2015l : IWalch fe Naabll2015h . Kimm fe Cenl 
( 203 ) argue that if more momentum is deposite d into a 
small er portion of the ISM (see, e.g., Ilffrig fe Hennebellel 
l2015h then this can lead to a larger amount of energy from 
supernovae being transferred to the ISM and galactic winds. 


The presence of a magnetic field does not appear to 
strongly affect the bulk properties of the resulting HII re¬ 
gion, even though it affects the structure of the cloud and 
HII region noticeably. If magnetic fields can alter the rate 
at which stars are formed by changing the structure of the 
densest parts of the cloud, then they would indeed alter the 
structure of HII regions by modifying the UV emission rate 
produced by the cluster embedded in the cloud. 

There are various pre-supernova stellar feedback pro¬ 
cesses that occur in star-forming clouds other than UV pho¬ 
toionisation. Radiation pressure from infrared and repro¬ 
cessed optical emission from stars can also aid cloud destruc¬ 
tion, although there is still some debate as to the effective - 
ness of these processes (e.g. iKrumholz fe Thompsonll2012i l. 
Radi ation pressure has recent ly been added to RAMSES- 
RT (iRosdahl fe Tevssiedl2015h and will be explored in fu¬ 
ture work. As well as radiative feedback, stars produce winds 
that can heat the gas arou nd them. Some theoretical wo rk 
has already been done by iGarcia-Segura fe Francol (Il996h 
in determining the effect these winds have on HII regions, 
particularly ultra-compact ones. For the majority of massive 
stars these winds are typically weak and provide a li mited 
amount of energy llDale et all 120141 : IGeen et al.l l2015l l , but 
more massive Wolf Rayet stars (> 20-30 Mg) do produce 
significant h ot bubbles in their circumstellar environment 
(|l)warkadasll2007i l . 

Another aspect not explored by these simulations is the 
role of metallicity, since we use only solar metallicity to cal¬ 
culate our cooling rates. Radiative cooling in primordial gas 
is particularly inefficient, and UV feedback at high redshift 
is expected to be even more efficient than in the simulations 
contained in this paper since the effects of metal cooling on 
the photoionised gas will be significantly reduced due to its 
low metallicity. Nonetheless, we expect the theoretical mod¬ 
els given in this paper to extend to HII regions with different 
metallicities. 
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10 CONCLUSIONS 


We present a new set of analytic models that de¬ 
scribe the evolution of HII regions in UV photo ionisa ¬ 
tion equilibrium ba s ed on t he arguments of|Spitze r| il978|); 

I Dyson fc Williams! lll98C fi : iFranco et al.l (Il990li : Matznerl 
1 2002ll ; Raga et^hl d2012l l . We focus on cases where the 
cloud’s density and velocity field corresponds to that ex¬ 
pected in a virialised cloud supported by turbulence. We 
present limits in which the emergence of HII regions from 
their host cloud is prevented by pressure forces, either by 
thermal pressure or by ram pressure from gas flows. In or¬ 
der to determine the validity of these ID models we compare 
them to a new suite of fully 3D simulations of UV photoion- 
sation feedback in turbulent, magnetised, self-gravitating 
clouds. We perform these simulations using RAMSES-RT, 
a Eulerian adaptive mesh refinement (AMR) radiation mag- 
netohydrodynamics code. 

The presence of turbulence is important for two com¬ 
peting reasons. On one hand, it provides a ram pressure term 
that resists the expansion of HII regions in cases where the 
source of UV photons is weak or the cloud is dense enough. 
We calculate an analytic estimate for the radius at which 
turbulence prevents the HII region from expanding further. 
If this radius is smaller than the radius of the cloud, the HII 
region remains trapped and does not destroy the cloud. On 
the other hand, turbulence provides support for the cloud 
and prevents it from collapsing. If a cloud is dominated by 
radial flows from infalling gas, the density increases dra¬ 
matically over time, requiring a much stronger source to 
destroy the cloud. Otherwise, the HII region collapses and 
is crushed by the cloud. Dense clumps orbiting the source 
position cause the HII region to “flicker” on and off as the 
dense gas efficiently absorbs the UV photons. 

The presence of a magnetic field limits fragmentation 
in the cloud, removing some channels of low-density gas 
through which the ionisation front can escape more easily, 
though the median radius of the ionisation front is similar 
in simulations with and without a magnetic field. 

We discuss briefly the expected role of UV photoionisa¬ 
tion in regulating the star formation efficiency (SFE) of the 
cloud. Since we do not directly model star formation in this 
work, we instead measure the mass in Jeans unstable gas, 
which we find decreases with increasing UV emission rate. 
We compare this to a simple model for the evaporation of 
dense clumps by UV photons, and find a reasonable agree¬ 
ment in most cases, though the model overpredicts the mass 
evaporated for the weakest photon emission rates. We will 
address self-consistent star formation and the amount of UV 
radiation from these stars in a future work. 
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APPENDIX A: FRONT EXPANSION IN A 
TURBULENT MEDIUM DERIVATION 

In this analysis we allow the velocity of the ionisation front 
and the shock to differ (in Section m we assume they are 
the same). We write the Rankine Hugoniot conditions (de¬ 
scribing the conditions across a shock interface) in the frame 
of the shock below, where v a is the shock velocity, and n c 
and v c are the post-shock number density and velocity (all 
other symbols as defined in Section [3] 


next (r, t) (v s - next (r, t)) = n c v c (Al) 


n ex t(r, t)((v a - next (r, t)) 2 + c 2 xt ) = n c (v 2 + c 2 x t) (A2) 


which can be combined to get: 


M e 2 xt = 

^ext 


v a - next (r, t) 


Cext 


(A3) 


hence, 


^ext 


v B - next (r, t)' 


(A4) 


Then from Raga et al. (2012), we can write the relation 
between the velocity of the ionisation front and the velocity 
of the shock: 


dry 

n s = — + n c 
d t 


and so 


giving 


dn 
Vs — — r~—h 


Vext 


d t Vs - next (r, t) 


(A5) 

(A6) 


%■ = (n s - next (r,t)) - Cext + n ext (r,t). (A7) 

Cl t Vs Vext [T , t) 

We now find an expression for (n s — n e xt(n, f)). As in 
Equation [3] photoionisation equilibrium gives 


2 3 o / / r,\2 2, 2 3 

riiVi * 3 / n e xt(n,0) r dr = n 0 r a . 
Jo 


(AS) 


Momentum conservation at the ionisation front gives: 

rue 2 = riext(r,t){vs - v ext {r,t)) 2 (A9) 

Using Equation IA81 this becomes: 

1 3/2 _ next (r, t) (v s - next {r,t)) 2 


and so 


Vs - next (r,t) 


no 


3/4 


no 


next (r, t ) 


1/2 


(A10) 


(All) 


Hence the equation for the velocity of the ionisation 
front is: 


1 dn 
d d t 


= F(r,t) - 


^ext 


next (r,t ) 


cf F(r, t) 


where 


F(r,t) = 


3/4 


( n 0 


n J \ next ( r , t) 


1/2 


(A12) 


(A13) 


APPENDIX B: CLOUD OUTFLOW MODEL 

We now discuss a model that takes into account the dif¬ 
ference in distance from the source to the cloud edge with 
angle, r c i oud . We sample r c ioud using rays along the vectors 
given in Appendix(D] We find that r c ioud has a roughly uni¬ 
form probability distribution function with limits at 3 and 
12 pc in the Fiducial cloud. 

We solve the radial evolution of the ionisation front up 
to v cloud using Equation^ with w = 0, i.e. a flat density field 
(since the spherically-averaged power law model no longer 
applies) with no set to the mean density of the cloud. We 
set /cloud to the time that the solution reaches r c i ou d. At 
this point, the ionisation front enters the diffuse external 
medium. We add a step function to this density profile where 
for r > r c i ou d, the density of the external medium is given 
by n d = 1 atom cm -3 . 

We modify Equation [2] to give 

ndn 2 =nicf. (Bl) 

We do not change Equation [3] since photoionisation equi¬ 
librium still holds, and none of these terms depends on the 
density of the external medium after t = 0. Solving as in 
Section m we find 

(t tcloud) — 

Aw f 1 + I V2 f ) 3/4 Ci(t ~ tcl ° Ud) 

y 4 \Tl d / \^cloud / ^cloud 

(B2) 

We solve this for r c ioud varying from 3 to 12 pc, giving a 
distribution of solutions for n. The median value has r c i ou d 
= 7.5 pc. We plot these solutions for each of the simulations 
in the Fiducial cloud in Figure [Bll 

The fit to this model is better than that of the power law 
for N48_B02 and N49_B02, which smooths over the sharp 
discontinuity at r c in cloud density. However, the model fails 
in the case of N48_B00, which has a density profile that is 
poorly fit by a flat field inside the cloud, where the “Power 
Law” model in Section ED is a better fit. The model also 
ignores clump motions and thus does not match N47_B02 
or N48_B02_C any better than the “Power Law” model. 
In addition, once the ionisation front has left the cloud its 
expansion no longer affects significantly the properties inside 
the cloud - for example, the momentum and ionised mass are 
well captured by the power law model. However, in cases 
where the expansion of the front once it leaves the cloud is 
important, this two-step outflow model is a useful tool in 
understanding the HII region. 



APPENDIX C: SHELL INSTABILITY 

In Section [8] we describe the loss of potentially star-forming 
mass from dense clumps as they are evaporated by the UV 
photons. However, as explained in Sections m and 16.21 the 
amount of mass in ionised gas is typically much smaller than 
the mass of neutral gas displaced by the ionisation front. 
Most of the mass displaced thus ends up in the dense, neu¬ 
tral shell around the ionisation front. Here we discuss the 
possibility that this shell becomes Jeans Unstable and thus 
able to fragment and form stars. 
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APPENDIX D: LINE-OF-SIGHT SAMPLING 

We sample lin es of sight on a sphere a ccording to the follow¬ 
ing algorithm feaff fe Kuiilaarslll997l f1. For N lines of sight 
with index i — {1, N}, the ith line of sight is defined as: 

, (2 i-N - 1\ _ , 

<pi = arccos I N _ ^ — 1 (Dl) 




3.6 \ 

sin 4>n\/N J 


(D2) 


where 9\ = 9 n = 0. The points on the surface of the sphere 
of unit radius n are then defined as 


fi 


sin (/>i cos 9i 
sin <j> i sin 9[ 
COS (f>i 


(D3) 


Figure Bl. Comparison of the “Outflow” model with the mean 
radius found in the simulation results as in Figure [6] Upper panel: 
varying photon flux. Lower panel: varying cloud compactness. 


lElmegreenl (|l994h states that a dense shell around a 
shock becomes unstable to fragmentation when 


nGpo 8 1 / 2 V r s heii 

3c o ffhell 


(Cl) 


where po is the initial density, Co is the sound speed in the 
shell, and T4heii and ?' a heii are the speed and radius of the 
shell respectively. We can assume these are n and n as be¬ 
fore. Differentiating Equation we get fi = tpri/t. For an 
initial density po = mn/X x 500 atoms/cm 3 with Co = 0.4 
km/s, we find that this criterion becomes 


r shell! > 100 pcMyr (C2) 


for the Fiducial cloud. This requires the shell to travel at 
least 25 pc in 4 Myr (roughly the lifetime of the most mas¬ 
sive star in our simulation) before it becomes unstable to 
fragmentation. We thus do not expect the shells around 
our ionisation fronts to become unstable in our simulations. 
However, on longer timescales it is possible that the shells 
would fragment and form stars. Another possibility is that 
this shell would encounter another molecular cloud, trigger¬ 
ing star formation via shock compression of this cloud. 


l 


see also http://people.sc.fsu.edu/~ jburkardt/f_src/sphere_grid/sphere_gri 



























